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ABSTRACT 


This  report  presents  survey  material  on  the  effects  of  the  reentry  environment  on  the  voltage 
breakdown  characteristic  of  antennas.  Enough  is  presently  known  about  these  effects  to 
indicate  their  relative  importance  and  to  allow  meaningful  experiments  to  be  designed  to 
diagnose  them.  Effects  of  convection  are  not  well  understood  because  of  conflicts  among 
various  theories  and  experimental  data.  High  gas  temperature  effects  are  not  adequately 
explained  and  insufficient  data  presently  exist  to  allow  extrapolation  with  confidence. 

Results  of  some  illustrative  calculations  of  flow  fields  about  a  typicL.  slightly  blunt-nosed 
body  are  presented,  including  an  assessment  of  the  vibrational  nonequilibrium  effect,  which 
is  found  to  be  important. 
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SECTION  1 


INTRODUCTION 

The  study  of  which  this  report  is  a  purt  is  concerned  with  the  effects  of  the  reentry  environ¬ 
ment  of  a  slender  hypersonic  vehicle  on  the  transmission  of  ECM  signals.  The  most  import¬ 
ant  aspect  of  this  study  is  the  effects  of  reentry  on  antenna  breakdown.  The  effort  during 
this  reporting  period,  as  well  as  the  last,  was  concentrated  on  a  survey  of  reentry  effects 
on  antenna  breakdown  and  included  vis  ts  to  the  facilities  of  some  of  the  workers  in  the  field. 
In  addition,  illustrative  calculations  of  matched  invisc id-boundary  layer  stagnation  line  flow 
field  calculations  have  been  performed  for  three  altitudes  of  reentry  of  a  0.5  inch  radius 
body,  and  the  effects  of  vibrational  nonequilibrium  in  the  shock  layer  at  250  kilofeet  altitude 
have  been  evaluated  by  a  superpositional  method. 

A  phenomenological  theory  of  reentry  effects  on  antenna  breakdown  is  surveyed  in  terms  of 
the  continuity  equation  and  the  net  ionization  frequency.  This  theory  is  applied  to  illustrative 
calculations  and  to  the  generation  of  constraints  on  future  experimental  parameters. 


J 
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SECTION  2 

FLOW  FIELD  ANALYSIS 


2.1  SHOCK  LAYER  STAGNATION  LINE  SOLUTIONS 

Profiles  of  the  flow  properties  along  the  stagnation  streamline  within  the  shock  layer  have 
been  obtained  from  matched  inviscid -boundary  layer  chemical  nonequilibrium  pure  air  solu¬ 
tions  for  altitudes  of  250,  200,  and  150  kiiofeet.  These  profiles  are  required  as  initial  con¬ 
ditions  for  the  nonequilibrium  boundary  layers  which  are  to  be  computed  during  the  last 
quarter  of  this  study;  they  are,  however,  of  considerable  interest  in  themselves  and  as  such 
are  presented  separately  from  the  boundary  layer  solutions. 

In  shock  layer  applications,  the  stagnation  line  matched  inviscid-boundary  layer  problem  is 
coupled,  i.  e. ,  the  inviscid  solution  must  be  used  to  provide  edge  conditions  for  the  boundary 
layer  solution,  which  in  turn  must  be  known  to  locate  the  bow  shock  as  influenced  by  bi  'indary 
layer  displacement.  The  matched  solution  must,  therefore,  be  obtained  by  iteration.  The 
procedure  used  herein  is  to  first  assume  the  boundary  layer  profiles  necessary  to  locate  the 
edge  at  which  the  inviscid  solution  is  terminated,  thereby  supplying  the  necessary  edge 
conditions.  The  boundary  layer  solution  with  these  edge  conditions  is  obtained  and  these 
profiles  are  then  used  to  relocate  the  edge.  The  process  is  repeated  until  the  change  in  the 
edge  is  less  than  a  specified  small  tolerance,  after  which  the  solutions  are  considered  matched. 

Before  proceeding  to  a  description  of  both  the  inviscid  and  boundary  layer  solutions,  it  is 
appropriate  to  point  out  that  both  solutions  have  identical  thermodynamic  and  nonequilibrium 
chemical  reaction  systems,  and  both  assume  the  condition  of  thermal  equilibrium. 

2.1.1  THE  INVISCID  SOLUTION 

In  the  inviscid  region  of  the  shock  layer,  the  solution  is  provided  by  the  streamtube  method 

of  McMenamin  and  O’Brien1  which  was  described  in  detail  in  the  First  Quarterly  Technical 

Report.  The  governing  variable  used  with  the  streamtube  method  is  the  streamwise  velocity 

2 

obtained  from  a  correlation  of  data  from  Burke,  Curtis  and  Boyer  ;  these  data  are  the  result 
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of  the  inverse  shock  nonequilibrium  inviscid  shock  layer  solution  of  Cornell  Aeronautical 
Laboratory,  'the  correlation  which  reproduces  the  data  of  Reference  2  to  within  +  8  %  is 
given  by: 

v  =  vg  (1  -  s/A)  (2.1.1.  -1) 

in  which  v  is  the  streamwise  velocity,  v  is  the  velocity  behind  the  shock  wave,  s  is  the 

s 

distance  from  the  shock  wave,  A  is  the  shock  detachment  distance,  and  the  exponent  is 
given  by: 

b  =  1  for  v  £  1670  ft/sec 

s 

b  =  0.  000596  v  for  v  >  1670  ft/sec 

s  s 

The  nonequilibrium  shock  detachment  distance  required  for  the  inviscid  shock  layer  solution 

is  obtained  in  the  following  manner.  By  examining  a  multitude  of  shock  detachment  distances 

from  equilibrium  air  and  perfect  gas  (lsys  5/3)  flow  field  solutions  for  spherical  noses, 
g 

Storer  found  that  all  cases  could  be  correlated  with  excellent  agreement  by  the  formula 


in  which  is  the  sphere  radius,  is  the  free  stream  density  and  p  is  the  density  behind 
a  l-jrmal  shock  wave  corresponding  to  either  the  equilibrium  air  or  perfect  gas  shock  solu¬ 
tion  as  the  case  may  be.  Since  the  gas  density  in  the  stagnation  region  for  either  equilibrium 
air  or  a  perfect  gas  is  nearly  constant,  it  might  be  expected  that  the  shock  detachment 
distance  could  be  essentially  specified  by  the  density  behind  the  norma]  shock.  This  cannot 
be  true  for  the  nonequilibrium  flow,  however,  because  the  density  varies  from  the  frozen 
value  (y  =  7/5)  at  the  shock  wave  to  the  equilibrium  value  at  the  stagnation  point.  Neverthe¬ 
less,  because  Storer's  formula  applies  equally  well  to  both  the  extremes  of  nonequilibrium 
flow,  it  might  be  expected  to  apply  also  in  the  nonequilibrium  case  if  the  density  behind  the 
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shock  were  replaced  by  a  suitable  average  density  along  the  stagnation  streamline.  Indeed, 
this  has  been  found  to  be  the  case  and  good  agreement  with  the  non  equilibrium  shock  detach¬ 
ment  data  of  Reference  2  has  been  obtained  whan  p  in  Storer's  formula  is  replaced  by 

s 


in  which  A  „  is  the  frozen  shock  detachment  distance,  and  p  is  the  density  from  the  non- 
FROZ 

equilibrium  streamtube  solution  which  is  obtained  using  the  velocity  correlation  with  A  = 

afroz* 

2.1.2  BOUNDARY  LAYER  SOLUTION 

In  the  boundary  layer  region  of  the  shock  layer  the  solution  is  provided  by  the  implicit  finite 

4 

difference  method  of  Blottner  .  At  the  stagnation  point  the  boundary  layer  equations  degen¬ 
erate  to  a  set  of  ordinary  nonlinear  differential  equations  with  two-point  boundary  conditions. 
After  a  transformation  to  the  Howarth-Dorodnitsyn  coordinates  system,  Blottner  solves  the 
linearized  set  of  governing  equations  by  an  implicit  finite  difference  scheme  which  requires 
the  successive  inversions  of  a  tridiagonal  matrix.  The  governing  equations  which  must  be 
solved  are  of  the  form 


+  a  —  +  a  H  +  a  =  0 
,2  Id  rj  2  3 

d77  ' 


(2. 1.2-1) 


in  which  rj  is  the  transverse  Howarth-Dorodnitsyn  coordinate  and  H  represents  the  longitudinal 
velocity  ratio  u/u  ,  the  temperature  ratio  T/T  ,  the  species  mass  fractions  C  ,  or  the  ele- 

6  ,  6  X 

ment  mass  fractions  CJ  depending  on  whether  the  momentum,  energy,  species  continuity,  or 

element  conservation  equation  is  being  considered.  Because  of  the  linearization  process  the 

coefficients  a  ,  a  ,  a  and  a  depend  on  the  unknown  variables  H.  Thus,  to  start  the 
1  Z  o  4 

solution,  initial  profiles  of  the  dependent  variables  must  be  provided  from  which  to  calculate 


the  coefficients.  The  set  of  governing  finite  difference  equations  forms  a  tridiagonal  matrix 
which  yields  the  solution  of  the  dependent  variables  after  inversion.  Based  on  these  values 
of  the  dependent  variables,  the  coefficients  are  recalculated  and  the  process  repeated  until 
the  change  in  the  dependent  variables  for  successive  iterations  is  less  than  a  specified 
small  tolerance. 

Calculations  for  a  0.5-inch-nose-radius  body  at  velocities  of  22. 1,  22. 0,  and  21. 95  kilofeet/ 
second  at  altitudes  of  250,  200,  and  150  kilofeet,  respectively,  have  been  done.  The  results 
of  the  matched  inviscid-boundaiy  layer  stagnation  line  solutions  are  presented  in  Figures 
2. 1.2-1  through  2.1.2-10,  where  successively,  ingroups  of  three,  the  inviscid  and  boundary 
layer  solutions  are  compared  at  altitudes  of  250,  200  and  150  kilofeet  for  the  gas  temperature, 
density,  electron  density,  and  pressure.  In  all  cases  the  boundary  layer  edge  was  arbitrarily 
defined  as  the  point  where  the  derivative  of  the  longitudinal  velocity  ratios  equaled  1/10,  i.  e. , 
at  the  point  where 

d  u/u 

— — 5  =  1/10  (2. 1.2-2) 

drj 

The  rather  large  mismatch  in  the  solutions  at  250  kilofeet  is  due  to  a  convergence  tolerance 
which  is  too  large,  an  effect  which  was  corrected  for  the  lower  altitudes.  In  all  three  cases 
the  boundary  layer  displacement  thickness  is  negative  due  to  the  cold  wall  condition  and,  with 
the  exception  of  the  highest  altitude,  it  is  so  small  as  to  be  unnoticeable  in  the  results.  Also, 
in  all  three  cases  the  catalytic  wall  condition  (thermodynamic  equilibrium)  is  imposed. 

It  should  be  noted  that  the  rapid  decrease  in  temperature  and  corresponding  increase  in 
density  for  the  inviscid  solution  near  the  wall  is  caused  by  the  equilibrium  condition  which 
is  intrinsically  required  by  reacting  inviscid  flows  at  a  stagnation  point. 
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NORMALIZED  DISTANCE  FROM  SHOCK  WAVE,  S/A 

Figure  2. 1.2-1.  Matched  Inviscid-Boundary  Layer  Solutions;  Gas  Temperature  vs  Normalized 

Distance  from  Shock  Wave  <250  kilofeet  Altitude) 
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GAS  TEMPERATURE  XlO-3,  T  (°K) 


NORMALIZED  DISTANCE  FROM  SHOCK  WAVE,  S/A 

Figure  2. 1.2-3.  Matched  Inviscid-Boundary  Layer  Solutions;  Gas  Temperature  vs  Normalized 

Distance  from  Shock  Wave  (150  kilofeet  Altitude) 
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GAS  DENSITY,  p  (SLUGS/FT3 


NORMALIZED  DISTANCE  FROM  SHOCK  WAVE,  S/A 

Figure  2. 1.2-4.  Matched  Inviscid-Boundary  Layer  Solutions;  Gas  Density  vs  Normalized 
Distance  from  Shock  Wave  (250  kilofeet  Altitude) 
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GAS  DENSITY,  p  (SLUGS /FT 


NORMALIZED  DISTANCE  FROM  SHOCK  WAVE,  S/A 

Figure  2. 1.2-6.  Matched  Inviscid-Boundary  Layer  Solutions;  Gas  Density  vs  Normalized 
Distance  from  Shock  Wave  (150  kilofeet  Altitude) 
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ELECTRON  NUMBER  DENSITY,  Ne  (cm“3) 


NORMALIZED  DISTANCE  FROM  SHOCK  WAVE,  S/A 

Figure  2.1  2-7.  Matched  Inviscid-Boundaiy  Layer  Solutions;  Electron  Density  vs  Normalized 

Distance  from  Shock  Wave  (250  kilofeet  Altitude) 
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ELECTRON  NUMBER  DENSITY,  Ne 


NORMALIZED  DISTANCE  FROM  SHOCK  WAVE,  S/& 

Figure  2. 1. 2-8.  Matched  Inviscid-Boundary  Layer  Solutions;  Electron  Density  vs  Normalized 

Distance  from  Shock  Wave  (200  kibfeet  Altitude) 
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NORMALIZED  DISTANCE  FROM  SHOCK  WAVE,  S/A 

Figure  2. 1.2-9.  Matched  Inviscid-Boundary  layer  Solutions;  Electron  Density  vs  Normalized 

Distance  from  Shock  Wave  (150  kilofeet  Altitude) 
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2 . 2  VIBRATIONA L  NONEQ UI  LIBRIUM 


The  relative  importance  of  vibrational  nonequilibrium  on  the  shock  layer  flow  properties 
has  been  investigated  by  applying  the  superpositioj  al  method  described  in  the  Second 
Quarterly  Technical  Report  to  the  inviscid  flow  field  at  250  kilofeet.  The  highest  altitude 
flow  field  was  selected  for  investigation  because,  since  vibrational  relaxation  times  vary 
inversely  with  pressure,  conditions  for  vibrational  nonequilibriurr.  are  most  favorable  at 
higher  altitudes.  Based  on  this  simplified,  order-of-magnitude  analysis  it  is  concluded  that 
vibrational  nonequilibrium  can  cause  substantial  changes  in  the  flow  properties,  particularly 
the  gas  temperature  and  density,  and  consequently  must  be  taken  into  account  to  accurately 
predict  the  shock  layer  properties  at  higher  altitudes.  The  analysis  and  the  results  which 
form  the  basis  for  this  conclusion  are  discussed  below. 


The  superpositional  method  used  in  this  investigation  consists  of  solving  the  uncoupled 
streamwise  vibrational  relaxation  equation  using  superposition  with  the  thermal  equilibrium 
chemical  nonequilibrium  inviscid  flow  field  to  provide  the  values  of  the  gas  pressure  and 
temperature  required  to  calculate  the  IocpI  relaxation  time  and  vibrational  equilibrium 

energy.  Only  the  primary  molecular  constituents  N  and  O  were  considered,  and  the 

*  “  5 

relaxation  times  for  each  of  these  species  were  obtained  from  Millikan  and  White  : 


-1 

1  (129T  '  -22.3) 
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-1 

1  (220T  1/0  -  24.  8) 

—  e 
P 


(2.2-1) 


(2.2-2) 


where  the  relaxation  times  Tq2  a°d  Tj^  ^re  in  seconds,  p  is  the  gas  pressure  in  atmospheres, 
and  T  is  the  gas  temperature  in  degrees  Kelvin.  The  harmonic  oscillator  model  is  used  to 
define  the  vibrational  temperature  corresponding  to  the  vibrational  energy  obtained  from  the 
solution  of  the  relaxation  equation  and,  therefore,  is  given  by 


(2.2-3) 
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=  9  /in  (1  + 
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R0  /e  ) 
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In  the  above  expression  and  are  the  vibrational  temperature  and  energy,  R  is  the  uni¬ 
versal  gas  constant  and  0^  is  the  characteristic  vibrational  temperature  which  has  a  value 
of  2239°K  for  oxygen  and  3395°K  for  nitrogen. 


Some  typical  results  of  the  analysis  are  shown  in  Figures  2.2-1  through  2.2-3  for  streamlines 
which  intersect  the  bow  shock  at  shock  angles  of  89,  70,  and  46  degrees,  respectively.  In 
these  figures  the  streamwise  variation  of  the  vibrational  temperatures  of  O  and  N2>  as  well 
as  the  gas  temperatures  obtained  from  solutions  for  thermal  equilibrium  chemical  nonequili¬ 
brium  and  for  perfect  gas  with  y  =  7/5  (undissociated  air  with  no  vibrational  excitation),  are 
given  as  functions  of  the  normalized  distance  along  the  body  surface.  In  all  cases  the  vibra¬ 
tional  temperatures  of  both  oxygen  and  nitrogen  are  characterized  by  an  initial  rapid  rise 
from  the  free  stream  value  at  the  shock.  However,  as  the  streamline  continues  to  expand 
around  the  body,  the  drop  in  the  gas  pressure  and  temperature  causes  a  rapid  increase  in 
the  local  relaxation  time  with  a  corresponding  rapid  leveling  off  ("freezing")  of  the  vibrational 
temperature.  If  at  this  point  the  vibrational  temperature  is  less  than  the  equilibrium  tempera¬ 
ture,  vibrational  temperature  will  continue  to  rise  gradually  until  the  two  are  equal,  which 
for  this  application,  particularly  for  nitrogen,  requires  a  considerable  distance  along  the 
body.  If  the  vibration  temperature  at  the  leveling  off  point  is  greater  than  the  equilibrium 
temperature,  the  vibration  temperature  will  decrease  gradually  until  the  two  are  equal.  It 
can  be  seen  that  this  latter  situation  is  the  case  for  the  89  degree  shock  angle  streamline, 
whereas  the  former  situation  applies  for  the  other  two  streamlines.  Also,  for  the  89  degree 
streamline  it  may  be  noticed  that  the  difference  between  the  no-vibrational-excitation  tempera¬ 
ture  and  the  thermal  equilibrium  temperature  back  on  the  conical  portion  of  the  body  is  about 
2500°K  or  about  45%  of  the  equilibrium  temperature.  In  this  case,  however,  this  temperature 
difference  should  not  be  interpreted  as  indicative  of  the  extent  of  vibrational  nonequilibrium 
effects  because  both  the  oxygen  and  nitrogen  vibrational  temperatures  are  above  the  thermal 
equilibrium  temperature.  As  a  result  of  this,  if,  in  the  solution  of  the  coupled  problem,  gas 
temperature  changes  from  thermal  equilibrium  do  not  significantly  affect  the  chemical  activity, 
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NORMALIZED  COORDINATE  ALONG  BODY  SURFACE,  Sfc/Rn 


Figure  2.2-1,  Vibrational  Nonequilibrium;  Streamwise  Variation  of  Temperature  vs 
Normalized  Coordinate  Along  Body  Surface  (Streamline  Shock  Angle  =  89  Degrees) 
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NORMALIZED  COORDINATE  ALONG  BODY  SURFACE,  Sb/R 

Figure  2,2-2.  Vibrational  Nonequilibrium;  Streamwise  Variation  of  Temperature  vs 
Normalized  Coordinate  Along  Body  Surface  (Streamline  Shock  Angle  =  70  Degrees) 
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NORMALIZED  COORDINATE  ALONG  BODY  SURFACE,  Sfc,/Rn 


Figure  2.2-3.  Vibrational  Nonequilibrium;  Streamwise  Variation  of  Temperature  vs 
Normalized  Coordinate  Along  Body  Surface  (Streamline  Shock  Angle  =  46  Degrees) 
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the  gas  temperature  from  the  coupled  solution  would  be  below  the  equilibrium  temperature 
shown  in  Figure  2. 1.2-1.  To  evaluate  the  extent  of  vibrational  nonequilibrium  effects  in 
this  case  would  require  a  comparison  with  the  gas  temperature  from  a  perfect  gas  solution 
with  y  =  9/7  (full  vibrational  excitation).  In  the  other  two  cases,  however,  where  the  vibra¬ 
tional  temperatures  are  below  the  thermal  equilibrium  temperature,  the  gas  temperature 
from  the  coupled  solution  would  be  between  the  y  =  7/5  curve  and  thermal  equilibrium  curve 
and,  consequently,  the  difference  between  these  two  curves  would  indicate  the  extent  of  the 
vibrational  nonequilibrium  effects. 
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SECTION  3 

REENTRY  EFFECTS  ON  ANTENNA  BREAKDOWN 


The  subject  of  this  section  is  the  effects  of  the  hypersonic  reentiy  environment  of  a  slender 
cone  body  on  voltage  breakdown  of  an  antenna  mounted  on  the  body.  The  discussion  surveys 
ideas  and  data  taken  from  written  reports  and  private  discussions  with  a  number  of  researchers 

g 

in  the  field.  As  indicated  in  a  previous  report  in  this  series  ,  the  various  slender  cone 
vehicle  flow  field  characteristics  which  can  affect  antenna  breakdown  may  be  listed  as  follows: 
convection,  gas  particle  number  density,  gas  temperature,  ionization,  chemical  additives, 
and  gradients  of  all  flow  field  quantities. 

The  prediction  of  the  effect  of  these  flow  field  characteristics  on  the  breakdown  field  of  a 
given  antenna  can,  in  principle,  be  accomplished  by  either  the  kinetic  theory  method  or  the 
phenomenological  method.  The  latter  method  is  more  instructive  in  a  general  survey  and 
also  it  can  fairly  readily  be  applied  to  reentry  phenomena  using  presently  available  theoretical 
and  experimental  data.  Hence  the  discussion  which  follows  treats  the  phenomenological  theory 
of  reentry  effects  on  antenna  breakdown.  Accordingly  the  treatment  will  follow  the  develop¬ 
ment  of  the  electron  continuity  equation,  which  describes  the  average  rates  of  electron 
number  density  gains  and  losses. 

Some  of  the  effects  of  reentry  enter  into  the  continuity  equation  directly  in  explicit  terms; 
others  are  implicit  in  various  terms  of  the  equation.  An  important  term  in  the  continuity 
equation  is  the  ionization  frequency  in  terms  of  its  dependence  on  the  microwave  field 
intensity,  which  allows  the  breakdown  field  to  be  related  to  the  solution  of  the  continuity 
equation.  Also,  the  antenna  near-field  distribution  has  an  important  impl'  it  effect  on  the 
solution  of  the  continuity  equation,  through  the  dependence  of  ionization  fre  ,aency  on  local 
field  intensity.  Local  gas  temperature  apparently  has  a  strong  effect  on  the  ionization 
frequency  at  the  temperatures  induced  by  hypersonic  reentry. 
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3.1  THE  CONTINUITY  EQUATION 

The  basic  starting  point  in  the  phenomenological  theory  of  antenna  breakdown  is  the  electron 
continuity  equation: 


V 


r 


dn 

dt 


(3.1-1) 


In  this  equation  the  electron  flux  vector  r  is  related  to  the  rate  of  change  of  electron  density 
n  with  time  t.  The  total  derivative  term  on  the  right-hand  side  is  meant  to  represent  only 
the  net  effects  of  sources  and  sinks,  such  as  ionization  and  attachment  rates  for  collisions 
of  electrons  with  neutral  particles.  If  only  electron-neutral  collisions  are  included  as 
electron  sources  or  sinks,  the  continuity  equation  can  be  written  in  the  following  form: 
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(3.1-2) 


where  v  is  the  average  net  ionization  frequency,  including  ionization  minus  attachment, 
n 

(We  omit  the  conventional  bar  over  the  symbol  for  average  frequencies  of  reactions,  since 
we  will  discuss  only  averages,  implicitly  assuming  that  the  effects  of  the  velocity  distribution 
can  be  properly  included  in  the  average  in  a  phenomenological  way. ) 

Before  one  can  proceed  to  the  next  step  in  the  analysis,  the  effect  of  convection  must  be 
considered.  The  fact  that  the  gas  which  is  being  subjected  to  the  high  near-fields  of  the 
antenna  is  moving  fairly  rapidly  past  the  antenna  means  that  either  the  electron  flux  must 
include  a  convective  term  or  the  frame  of  reference  for  time  variations  must  move  with  the 
flowing  gas.  These  two  alternative  views  of  convection  give  rise  to  the  two  prevalent 
classes  of  convection  theories. 


3.1.1  FINITE  TIME  THEORY 

In  the  frame  of  reference  of  an  observer  moving  with  the  gas  velocity  past  the  antenna,  the 
local  electric  field  intensity  appears  to  vary  with  time  in  the  same  way  as  the  variation  of 


antenna  near-field  along  a  streamline  of  the  flow.  In  the  very  simple  situation  of  uniform 
flow  with  velocity  V  through  a  region  of  uniform  field  of  length  L,  the  gas  will  be  affected 
by  the  field  for  a  length  of  time  equal  to  L/V.  In  this  case  the  analysis  can  be  performed 
in  exactly  the  same  way  as  pulsed  signal  breakdown  of  a  nonflowing  gas,  where  the  pulse 
duration  is  taken  as  L/V.  (Of  course,  for  pulsed  signals  in  such  a  situation  the  effective 
pulse  duration  is  either  L/V  or  the  actual  pulse  duration,  whichever  is  shorter. )  This  is 

7 

Kelly  and  Margenau's  method  .  The  preceding  simple  finite  time  result  is  not  strictly 
rigorous  even  under  the  simple  uniform  field,  uniform  flow  conditions  postulated,  since  the 
diffusion  losses  in  the  presence  of  the  flow  are  not  the  same  as  those  in  the  absence  of  flow. 
For  example,  in  a  uniform  field  region  bounded  by  parallel  plates,  the  electron  density 
distribution  is  symmetrical  about  the  midplane  (the  lowest  order  mode  gives  the  well  known 
result  for  the  diffusion  length  L/rr).  But  as  shown  in  Figure  3. 1. 1-1,  in  the  presence  of  gas 
flow  between  the  plates  the  electron  density  distribution  must  be  unsymmetrical,  with  the 
peak  shifted  toward  the  downstream  side,  since  the  gas  on  the  downstream  side  has  been 
under  the  influence  of  the  field  for  the  longest  period  of  time.  Note  that  this  case  is  discussed 
only  to  illustrate,  in  the  most  fundamental  way  possible,  how  convection  affects  diffusion,  so 
that  the  finite  time  result  does  not  give  the  complete  answer.  We  do  not  wish  to  imply  that 
the  situation  depicted  in  Figure  3. 1. 1-1  is  representative  of  the  reentry  problem. 

The  above  objection  to  the  finite  time  theory  of  convection  i3  fundamentally  related  to  the 
problem  of  coupling  between  convection  and  diffusion.  Thu  objection  alone  may  not  be  very 
serious  except  in  terms  of  the  Cottingham  theory  (see  below),  but  it  will  be  true  that  the 
reentry  flow  field  environment  makes  it  particularly  difficult  to  follow  the  breakdown  process 
in  the  frame  of  reference  of  the  flowing  gas.  Velocity  gradients,  or  shear  effects,  will 
unduly  complicate  the  problem  mathematically.  Also,  the  antenna  nonuniform  near-field 
distribution,  even  in  the  presence  of  uniform  flow,  would  present  a  difficult  problem,  since 
it  would  be  analogous  to  the  problem  of  breakdown  of  a  nonrectangular  pulse. 
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ELECTRON  DENSITY,  n/npeak 


Figure  3. 1. 1-1.  Effect  of  Gas  Flow  on  the  Electron  Density  Distribution  in 
a  Uniform  Field  Region,  with  Flow  Normal  to  the  Walls 


3.1.2  COUPLED  CONVECTION-DIFFUSION  THEORY 

The  class  of  convection  theories  which  couple  the  convection  and  diffusion  phenomena  is 
based  on  the  use  of  the  following  equation  for  the  electron  flux  vector: 

r  =  nV-V  (Dn)  (3. 1.2-1) 

where  V  is  velocity  and  D  is  the  diffusion  coefficient.  When  this  ti.  .session  is  sub¬ 
stituted  in  Equation  3. 1.2-2,  the  continuity  equation  becomes 

an  2 

?>t  ~  Vnn  +  V  (Dn)"v,(nV)  (3. 1.2-2) 

_  8 

Covert  has  suggested  a  change  of  variable  which  eliminates  the  term  in  the  divergence 
from  the  continuity  equation.  This  change  of  variable  may  be  generalized  to  allow  for 
spatial  variation  of  the  diffusion  coefficient. 

Let 
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>b  exp  ( 
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•  dr) 


(3. 1.2-3) 


where  r  is  the  radius  vector  from  an  arbitrary  origin.  The  exponent  in  this  expression 
is  proportional  to  the  integral  along  a  streamline  of  the  ratio  of  the  diffusion  time  to  the 
convection  time.  The  resulting  form  of  the  continuity  equation  is 
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>/> 


1  9n 
Dn  dt 


(3.  1,2-4) 


If  separation  of  variables  in  space  and  time  can  be  accomplished,  this  equation  represents 
a  homogeneous  differential  equation  with,  in  general,  a  nonconstant  coefucisnt. 
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Compare  the  form  of  Equation  3.1.  2-4  with  the  equation  which  would  result  if  convection 
were  absent. 
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V 


(Dn)  + 


1  ^n 
Dn  dt 


Dn 
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(3. 1.2-5) 


This  equation  is  of  exactly  the  same  form,  implying  that  the  net  ionization  frequency  with 
convection  is  given  by 


(3. 1.2-6) 


where  the  subscript  o  represents  the  absence  of  convection.  This  result  actually  can  be 
true  only  at  one  point  in  the  case  of  nonuniform  electric  field,  but  then  it  can  be  interpreted 
in  terms  of  a  breakdown  criterion  at  the  peak  field  point.  With  the  exception  of  the  more 
general  term  permitting  compressible  flow.  Equation  3- 1.2-6  is  identical  with  the  results 
of  Romig  and  Fante  ,  which  were  derived  without  the  aid  of  Covert's  transformation,  but 
with  the  use  of  boundary  conditions  and  separation  of  variables.  Thus  the  transformation 
allows  the  coupling  between  diffusion  and  convection  to  be  displayed  directly  in  the  form  of 

the  ratio  V/(2D)  without  the  necessity  of  invoking  boundary  condition  or  separation  of 

» 

variables , 


The  example  cited  above  of  the  simple  one-dimensional  situation  of  uniform  gas  flow 
between  two  parallel  plates  bounding  a  uniform  field  region  offers  a  simple  illustration  of 
the  result  of  the  coupled  convection-diffusion  theory.  The  lowest  order  solution  of  Equation 
3. 1.2-4,  for  the  condition  of  vanishing  electron  density  at  the  boundaries  and  the  steady 
state  condition  of  vanishing  di/j/dt,  is  the  following: 

0  =  >/<  sin  —  (3. 1.2-7) 

1  L 
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(3. 1.2-6) 


where  0  is  the  maximum  value  of  ijj  and  x  is  the  distance  between  the  plates.  Thus  the 
electron  density  distribution  is 


n  = 


n^  exp 


(X-X1>  2D 
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(3. 1.2-9) 


sin- 


where  n^  is  the  peak  electron  density,  occurring  at  the  point  x^. 
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tan 
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(3.1.2-10) 


This  result  has  been  plotted  in  Figure  3. 1. 1-1  for  the  case  of  x^  =  3L/4.  The  value  of 

V  for  this  case  is  2  7r  D/L,  which  would  be  on  the  order  of  20,  000  feet/second  for  air  near 

the  breakdown  minimum  using  tf  a  frt  diffusion  coefficient.  Also,  the  value  of  \>  is  twice 

n 

its  value  without  flow  in  this  exam^.e. 

The  above  simple  calculation,  besides  illustrating  how  the  theoiy  is  applied,  shows  the 
magnitude  of  the  results  which  it  predicts.  It  will  be  seen  below  that  the  breakdown  field 
varies  as  a  relatively  low  power  of  the  ionization  frequency,  so  that  the  convection  of  air 
at  20, 000  feet/second  is  predicted  to  increase  the  minimum  breakdown  power  by  only  about 
1  db.  Of  course,  a  greater  effect  would  be  observed  in  the  presence  of  ionization,  since  the 
ambipolar  diffusion  coefficient  is  much  less  than  the  free  diffusion  coefficient,  but  it  will 
be  seen  that  the  breakdown  criterion  then  becomes  an  important  consideration  which  modifies 
the  result. 
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There  are  a  few  fine  points  which  have  been  left  out  of  the  above  discussion,  since  it  was 

meant  only  as  a  very  simple  example.  Equation  3. 1. 2-8  follows  directly  as  a  condition 

for  the  existence  of  a  solution;  but  Equation  3. 1.2-9  contains  an  arbitrary  constant,  the 

10 

peak  electron  density,  implying  that  the  solution  is  not  unique.  The  solution  of  Fante 
for  the  problem  of  breakdown  between  parallel  plates  with  flow  parallel  to  the  plates  is 
just  the  opposite.  There  the  electron  density  distribution  is  a  unique  function  of  the  net 
ionization  frequency;  and  the  value  of  given  by  Equation  3. 1.2-8  is  a  threshold  above 
which  the  solution  fails  to  exist.  Of  course,  the  one-dimensional  flow  problem  solved 
above  is  physically  impossible  to  realize  without  having  the  boundaries  be  perfectly  porous 
to  all  particles  but  electrons.  However,  it  may  approximate  the  situation  in  which  the 
afterglow  decay  beyond  a  uniform  field  region  is  extremely  rapid. 

W.  C.  Taylor  (Stanford  Research  Institute,  private  communication)  has  pointed  out  that 
Fante' s  result  applies  only  if  the  breakdown  criterion  is  that  the  electron  density  becomes 
infinite.  He  finds  that  the  peak  electron  density  at  lower  values  of  Vr  is  actually  increased 
by  orders  of  magnitude  over  the  initial  value, and  that  the  finite  time  theory  is  a  good  approxi¬ 
mation  to  the  exact  results  if  the  breakdown  criterion  is  relaxed  to  require  only  a  significant 
rise  from  the  initial  electron  density.  This  result  makes  the  convection  effect  even  weaker 
than  it  has  been  represented  above.  In  other  words,  the  coupled  convection-diffusion  theory 
represents  an  upper  limit  for  the  net  ionization  frequency  at  breakdown, and  the  finite  time 
theory  is  a  good  approximation  if  a  finite  electron  density  is  used  as  a  breakdown  criterion. 
Private  communications  with  various  investigators  reveal  that  this  viewpoint  apparently 
agrees  with  shock  tube  experiment  results  for  single  pulse  breakdown.  These  latter  considera¬ 
tions  probably  apply  only  to  the  case  of  significant  flow  field  ionization,  which  causes  ambi- 
polar  diffusion  and  the  modification  of  the  breakdown  criterion.  But  convection  has  a  minor 
effect  in  the  presence  of  free  diffusion. 

3. 1.2.1  Space  Change  Coupled  Convection 

The  preceding  discussion  might  lead  to  the  conclusion  that  the  convection  effect  is  relatively 
unimportant,  except  for  a  theory  of  W.  B.  Cottingham  (Purdue  University,  private 
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communication).  This  theory  takes  account  of  the  coupling  between  space  charge  field  and 
convection  and  diffusion  of  both  ions  and  electrons.  It  is  very  similar  in  form  to  the  simple 
coupled  convection-diffusion  theory  discussed  above,  but  it  predicts  a  much  different  result 
for  the  effect  of  convection  on  breakdown  field.  Cottingham  has  verified  this  theory  in 
experiments  on  multiple  pulse  breakdown  in  pipe  flow. 

The  derivation  starts  with  expressions  for  the  electron  and  ion  flux  which  resemble  those 
of  Allis  and  Rose11. 


T  =  nV  -  v  (D_n)  -  (iEn  (3.1. 2.1-1) 

I\  =  n  V-v  (D  n J  +  u  En^  (3. 1.2. 1-2) 

'  T  +  +  I  “T 


where  D_  is  the  electron  free  diffusion  coefficient,  p  is  the  electron  mobility,  E  is  the 
space  charge  field,  and  the  subscript  +  refers  to  the  positive  ions.  In  the  steady  state  the 
electron  flux  is  equal  to  the  ion  flux  and  the  space  charge  field  may  be  eliminated  between 
these  equations: 


r  = 


(p+p  +  )nn+V-  P+n+V(D_n)  -  pnV(D  +  n+) 
pn  +  p  +  n  + 


(3. 1.2. 1-3) 


According  to  Allis  and  Rose,  the  free  diffusion  of  electrons  at  low  electron  densities  is 
equivalent  to  the  condition 


D  n  =  D  n 
+  + 


(3. 1.2. 1-4) 


\ 


\ 

\ 


\ 
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mmM 


which  gives 


r  = 


^  +  M, 


«D++  „+D_  D-nV  ‘  7  <D-"> 


-V- 


/  nV-V  (D_n) 


(3. 1.2. 1-5) 


where  is  the  ambipolar  diffusion  coefficient.  Also,  according  to  Allis  and  Rose,  the 
ambipolar  diffusion  of  electrons  at  high  electron  densities  corresponds  to  equal  electron 
and  ion  densities,  giving 


r  «  nV  -  — - —  v  (D  n) 

p  +  p+  - 


_JL_ 
H  +  M+ 


7  (D+n) 


(3. 1.2. 1-6) 


Under  the  reasonable  assumption  that  the  ratio  of  electron  and  ion  mobilities  is  independent 
of  position,  this  gives 


r  =  nV  -  V  (D  n) 
a 


(3. 1.2. 1-7) 


It  follows  from  Equations  3. 1. 2. 1-5  and  3. 1. 2. 1-7  that  the  electron  flux  vector  is  given  by 
the  following  equation  at  the  limits  of  either  free  or  ambipolar  diffusion. 


r 


nV  -  V(Dgn) 


(3. 1.2. 1-8) 


where  Dg  is  the  transitional  electron  diffusion  coefficient.  This  is  Cottingham's  equation. 
Under  the  "constant  ratio  approximation"  (see  below)  this  equation  would  apply  to  the 
transition  from  free  to  ambipolar  diffusion  as  well  as  to  just  either  limit,  as  derived  above. 
The  true  form  of  this  equation  can  be  determined  only  by  simultaneous  solution  of 
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Equations  3. 1. 2. 1-1  and  3,1.2. 1-2  with  Poisson's  equation,  as  done  by  Allis  and  Rose 

12 

for  the  case  of  no  convection;  but,  as  pointed  out  by  Allis  ,  an  equation  like  3. 1.  2. 1-8 
is  approximately  correct  when  the  transitional  diffusion  coefficient  is  related  to  the  electron 
density  by  the  proper  criterion. 


It  is  immediately  apparent  that  the  last  equation  is  of  the  same  form  as  Equation  3. 1.  2-1. 
Thus  the  space  charge  coupled  convection  theory  can  be  viewed  as  a  generalization  of  the 
simple  diffusion  coupled  convection  theory  given  above.  Then  the  generalized  continuity 
equation  is 


(3. 1.2. 1-9) 


tj)  =  Dgn  exp 


(3.1.2.1-10) 


where  V  is  to  be  interpreted  as  the  flow  velocity  v  in  the  ordinary  coupled  diffusion- 

convection  theory  and  as  D  v/D  in  the  coupled  space  charge  theory. 

s  a. 

A  simple  physical  interpretation  of  the  latter  theory  can  be  made  as  follows.  The  convection 
of  ions  in  the  background  of  flowing  gas  creates  a  space  charge  force  on  the  electrons  which, 
because  of  the  greater  mobility  of  the  electrons,  makes  convection  more  efficient  by  the 
ratio  of  the  (transitional)  diffusion  coefficient  to  that  for  ambipolar  diffusion.  In  fact  if  it 
is  noted  that  the  gas  convection  acts  directly  only  on  the  ions,  because  the  electron  energy  . 
far  exceeds  the  kinetic  energy  of  convective  motion,  then  the  convective  term  in  Equation 
3. 1. 2. 1-1  should  be  dropped.  The  first  term  *n  tht  rumerator  of  Equation  3. 1. 2. 1-3  would 
then  involve  only  q  and  not  also  jt+,  but  this  involves  no  difference  in  the  results,  since  ju  is 
much  greater  than  u  and  is  negligibly  different  from  their  sum. 

"P 
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Now  the  predicted  result  of  convection  is  not  negligible,  especially  when  the  the  diffusion 
is  free,  since  Dg/Da  can  be  as  large  as  100.  This  means  that  the  result  shown  in  Figure 

3. 1.1- 1  can  correspond  to  a  velocity  on  the  order  of  only  200  feet/second.  Also,  Equation 

3.1. 2- 8  in  the  space  charge  coupled  theory  takes  the  form 


l—\2  + 

/  *  I' 

K) 

(3.1.2.1-11) 

And  Equation  3. 1.2-6  is 


_V_\ 

2D  )  (3.1.2.1—12) 

a  / 

Therefore,  a  velocity  of  20,  000  feet/second  could  increase  the  ionization  frequency  by  a 
factor  of  10  .  But  if  this  much  increase  in  ionization  frequency  were  to  occur,  it  would 
probably  create  ambipolar  diffusion,  in  which  case  Equation  3.1.2.1-12  would  be,  neglect¬ 
ing  the  divergence, 

D  2 

_  _ a  V 

Vn  vn  D  4D  (3.1.2.1-13) 

°  -  a 


where  is  the  ionization  frequency  without  flow  under  free  diffusion  conditions.  This 
predicts  that  a  flow  velocity  of  20, 000  feet/second  would  increase  the  ionization  frequency 

by  a  factor  on  the  order  of  100,  which  corresponds  to  about  4  db  increase  in  the  minimum 
breakdown  power  of  air. 

3. 1.3  SEPARATION  OF  VARIABLES 

Up  to  this  point  the  time  derivative  of  the  electron  density  has  been  retained  in  the  expression 
for  the  coefficient  of  the  unknown  function  in  the  continuity  equation.  The  discussion  in  the 
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previous  section  was  based  on  the  steady  state  assumption  of  vanishing  time  derivative, 
which  corresponds  to  continuous  wave  (CW)  breakdown.  In  order  to  obtain  a  solution  for 
pulsed  signal  breakdown  it  is  necessary  to  separate  variables  with  respect  to  space 
and  time. 


It  is  customary,  though  not  always  explicitly  stated,  to  assume  that  the  quantities  v  and 
Dg  Equation  3. 1.2. 1-9  are  independent  of  time  and,  at  most,  functions  only  of  position. 
This  assumption  is  implicitly  equivalent  to  assuming  that  the  electron  energies  relax  to 
their  steady  state  distribution  much  more  rapidly  than  the  electron  density  distribution  does. 
The  t’me  independence  of  also  implies  that  the  electron  density  does  not  vary  sufficiently 
widely  to  change  the  value  of  the  transitional  diffusion  coefficient  over  an  appreciable  fraction 
of  the  time  for  onset  of  breakdown.  Of  course,  V  is  an  external  parameter  determined  by 
the  flow  field  properties  and  it  is,  at  most,  a  function  of  position.  Equation  3. 1.  2. 1-9 
separates  as  follows: 


D  -2JS 
s  j b 


—  d  n 
n  St 


(3. 1.3-1) 


where  y  is  a  separation  constant. 


n 


n  e 
o 


y  t 


(3. 1.3-2) 


where  nQ  is  the  initial  electron  density.  Also  for  a  pulsed  signal 


(3. 1.3-3) 


where  is  the  electron  density  at  breakdown  and  r  is  the  pulse  time.  The  continuity 
equation  in  spatial  coordinates  is  now 
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(3. 1.3-4) 


Since  y  vanishes  in  CW  breakdown,  the  above  result  shows  that  the  breakdown  criterion 
for  the  ionization  frequency  for  pulsed  signals  is: 


v  =  v 

pulse  nCW 


1 

T 


Jin 


'B 


n 


(3. 1.3-5) 


which  is  the  result  found  by  Gould  and  Roberts 
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3.1.4  EFFECTS  OF  THE  REENTRY  ENVIRONMENT 

The  continuity  equation  generally  leads  naturally  to  a  breakdown  criterion  involving  the 
bracketed  terms  comprising  the  coefficient  of  $in  Equation  3. 1.3-4  in  addition  to  a  term 
involving  a  characteristic  length  determined  by  the  geometry  of  the  particular  problem. 

The  effects  of  the  reentry  environment  appear  implicitly  in  each  of  these  terms,  besides  the 
already  discussed  explicit  appearance  of  the  convective  velocity.  Now  consider  each  reentry 
phenomenon  in  terms  of  its  effect  on  these  terms. 

3.1. 4.1  Gas  Density 

The  classical  theory  of  microwave  breakdown  expresses  breakdown  criteria  in  terms  of 
proper  variable  such  as  v^/p  and  Dp,  where  p  is  the  pressure  in  torrs.  These  quantities 
generally  normalize  in  approximately  this  way  because  of  the  dependence  of  collision  fre¬ 
quencies  on  number  density  of  gas  particles.  In  voltage  breakdown  experiments  under 
standard  laboratory  conditions,  the  practice  is  to  express  the  state  of  the  gas  being  tested 
in  terms  of  pressure,  since  this  is  a  directly  measurable  quantity.  However,  pressures  in 
the  reentry  induced  shock  layer  are  often  orders  of  magnitude  higher  than  in  the  free  stream. 

Also,  because  of  increased  gas  temperature,  shock  layer  density  is  not  generally  the  same 

14 

as  In  the  free  stream.  Whitmer  and  MacDonald  have  noted  that  the  physically  important 
parameter  for  breakdown  is  the  gas  density  rather  than  pressure,  since  it  is  the  number  of 
gas  particles  per  unit  volume  which  determines  the  rates  for  ail  the  processes  which  control 
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breakdown.  Thus  the  use  of  either  the  free  stream  pressure  or  the  shock  layer  pressure 
in  predicting  breakdown  at  a  given  altitude  would  give  erroneous  results.  The  distribution 
of  gas  density  in  the  high  field  intensity  region  is  the  correct  determining  factor. 

In  order  to  give  a  quantity  which  is  commensurable  with  the  classically  used  parameter 
of  pressure,  most  authors  define  a  reduced  pressure  p*.  Reilly  defines  the  reduced 
pressure  in  torrs  in  terms  of  the  gas  number  density  N  in  particles  per  cubic  centimeter: 

p*  =  2.9x10~17N  (3. 1.4. 1-1) 

16 

Light  and  Taylor  define  the  reduced  pressure  in  an  equivalent  way  in  terms  of  tempera¬ 
ture  T  in  °K: 

p*  =  273  p/T  (3. 1.4. 1-2) 

These  corrections  have  the  effect  of  shifting  the  curve  of  breakdown  field  versus  altitude 

about  30, 000  feet  higher  in  altitude  at  velocities  of  the  order  of  22,  000  feet/second,  accord- 

14 

ing  to  Whitmer  and  MacDonald  . 

1 7 

Gas  density  effects  have  been  observed  in  reentry  flight  experiments.  Bisbing  showed 
that  the  observed  period  of  breakdown  of  a  slot  antenna  on  a  conical  vehicle  agreed  with 
laboratory  data  in  cold  air  when  the  air  density  corresponding  to  the  inviscid  flow  region 
was  used  as  the  independent  variable.  Nanevicz  also  showed  that  when  the  vhf  breakdown 
data  from  Nike-Cajun  flights  is  corrected  for  local  density,  it  agrees  roughly  with  laboratory 
measurements  in  cold  air. 

The  gas  density  affects  all  the  terms  of  the  continuity  equation  in  their  proper  variables 
form  in  terms  of  the  reduced  pressure  p*,  but  the  remaining  reentry  environmental  factors 
affect  only  one  or  two  terms  each.  Therefore,  the  remainder  of  the  discussion  of  the  effects 
of  reentry  on  the  continuity  equation  treats  each  of  these  terms  in  turn. 
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3.1.4.  2  The  Diffusion  Coefficient 

It  is  conceivable  that  the  diffusion  coefficient  will  be  affee  ed  by  flow  field  ionization,  gas 

temperature,  and  chemical  composition.  The  ionization  effect  on  the  diffusion  coefficient 

is  well  known  and  has  been  discussed  by  many  authors.  As  indicated  above,  it  is  necessary 

to  recognize  this  fact  by  writing  the  diffusion  coefficient  as  D  ,  denoting  transitional  dif- 

s 

fusion.  The  magnitude  of  the  transitional  diffusion  coefficient  can  be  calculated  only  by 
solving  the  boundary  value  problem  constituting  the  separate  electron  and  ion  continuity 
equations,  with  space  charge,  and  Poisson's  equation.  The  result  depends  on  the  geometry 
of  the  boundaries  as  well  as  on  the  ionization  level. 


Allis  and  Rose  have  performed  calculations  for  hydrogen  and  found  that  the  transition 

from  free  to  ambipolar  diffusion  takes  place  over  five  orders  of  magnitude  of  electron 

density.  Moreover,  their  results  indicate  that  the  constant  ratio  approximation  is  nearly 

correct  if  the  electron  density  used  in  this  approximation  is  decreased  from  its  actual  value 

by  a  constant  factor.  From  Figure  4  of  their  paper  it  appears  that  this  factor  should  be 
12 

about  20,  but  Allis  indicates  a  factor  of  five  in  another  calculation.  A  factor  of  ten  may 

be  reasonable.  The  constant  ratio  approximation  consists  of  assuming  that  the  ratio  of 

12 

electron  and  ion  densities  is  constant  throughout  the  discharge.  Allis  has  derived  an 
expression  for  the  transitional  diffusion  coefficient  which  follows: 


D 

s 

D  D 
-  a 


(3.1.4. 2-1) 


where  A  is  the  diffusion  length  (see  below),  £  ^  is  the  Debye  length,  p  is  mobility,  n  is 

particle  density  (+  refers  to  ions),  D  is  the  electron  free  diffusion  coefficient  and  D  is 

a 

the  ambipolar  diffusion  coefficient: 
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(3. 1.4. 2-2) 


/Li  O  +  /LiD 

D  =  _L_z - + 

a  /u+  +  ,u 


l 


e  D 
o  - 

ne/.' 


(3. 1.4. 2-3) 


where  e  is  the  permittivity  of  free  space  ande  is  the  electronic  charge „  T~  oractice,  the 
ratio  ^i+  n+  / (/un)  is  generally  small  compared  with  unity,  except  possibly  when  the  Debye 
length  is  very  large,  in  which  case  the  ratio  is  not  important.  Thus  a  good  approximation 
which  depends  only  on  the  electron  density  is  the  following  equation. 


D 

_ s_ 

D  D 


-  a 


,  2  2 

£d  +/ 


Da  V  +  D-  A" 


(3. 1.4. 2-4) 


In  the  limits  of  low  and  high  electron  density  the  value  of  Dg  approaches  D  and  D  respectively. 
In  order  to  allow  the  above  constant  ratio  approximation  formula  to  reproduce  approximately 
the  actual  transition  situation,  apply  the  factor  of  ten  recommended  above. 


D 

_ s_ 

D  D 


a 


10  XD2  +  A2 

10  D  ^2  +  D  A2 
a  D 


(3. 1.4. 2-5) 


Thus  the  transition  from  free  to  ambipolar  diffusion  in  this  approximation  occurs  when  the 
Debye  length  is  about  one  third  of  the  diffusion  l3ngth. 

Two  simplified  approximations  to  the  ambipolar  diffusion  coefficient  are  possible,  considering 
the  fact  that  n  is  much  greater  than  n+.  In  an  active  discharge,  D  is  much  greater  than 
D+p,  so  that  the  ambipolar  diffusion  coefficient  is  given  approximately  by  the  following  equation: 
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But  in  the  case  of  an  isothermal  plasma,  D  n  is  equal  to  D  n  : 

—  +  4* 


2,+ 

D  =  -  D  =  2D 

a  + 


'3. 1.4. 2-6) 


(3.1.4. 2-7) 


It  is  important  to  note  that  the  latter  result  applies  only  to  the  afterglow  plasma,  for  example 
the  decaying  plasma  which  exists  between  pulses  in  multiple  pulse  breakdown.  Also,  there 
is  not  actually  a  factor  of  two  difference  between  these  two  results,  because  D  / \i  is  much 
greater  in  the  active  discharge  than  in  the  isothermal  plasma,  while  n+  is  practically  constant. 

It  is  apparent  that  one  needs  values  for  the  parameters  D  and  ju+  / \i  to  apply  these  equations 
and  also  one  needs  a  value  of  D  / \l  in  order  to  calculate  the  Debye  length  for  Equation 
3. 1.4.  2-5.  The  latter  problem  is  accurately  enough  handled  for  the  purposes  of  the  transi¬ 
tional  diffusion  in  the  constant  ratio  approximation  by  casting  Equation  3. 1.4. 2-3  in  the 
following  form: 


e  kT 
o  e 

2 

ne 


(3. 1.4.2 -8) 


where  k  is  Boltzmann's  constant  and  T  is  the  electron  temperature.  In  the  isothermal, 

afterglow  plasm?  T  is  equal  to  the  gas  temperature;  and  in  an  active  discharge  in  air  kT 
e  e 

is  on  the  order  of  a  few  electron  volts.  At  breakdown  in  air,  a  commonly  used  value  of  D 

7 

is  given  by  Kelly  and  Margenau  : 
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But  MacDonald  prefers  to  use  a  value  which  depends  on  the  microwave  field  intensity  E. 

10  ‘D  p  =  29+0.  9  E  /p  cm^  torr  sec  1  (3.1.4.2-10) 

e 

where  E^  is  the  effective  field  (see  below). 

Various  values  of  n/n  (or  D  /D  )  for  breakdown  in  air  have  appeared  in  the  literature, 

+  “  7a  14 

including  Kelly  and  Margenau's  value  of  100  and  Whitmer  and  MacDonald’s  value  of  40. 

But  these  are  not  actually  inconsistent,  since  the  latter  authois  applied  Equation  3. 1.4.  2-7 

to  the  data  on  mobilities  for  nonisothermal  plasmas.  Tims  the  value  of  40  is  actually  80 

and  it  seems  reasonable  to  assume  that  for  breakdown  of  air: 

D 

~  2=100  (3.1.4.2-11) 

D  u 

a  + 

The  discussion  of  the  diffusion  coefficient  so  far  has  concentrated  on  the  effects  of  flow 
field  ionization.  It  was  mentioned  above  that  gas  temperature  and  composition  can  conceivably 
affect  diffusion.  The  gas  composition  effects  are  fairly  obvious,  it  being  apparent  that  mole¬ 
cular  dissociation  of  the  air  and  the  formation  of  products  of  ablation  in  the  flow  field  sur¬ 
rounding  a  reentry  vehicle  would  require  one  to  obtain  basic  data  on  the  mobilities  of  electrons 

and  ions  in  the  product  gases,  some  of  which  are  undoubtedly  available  in  the  literature.  The 

20 

gas  temperature  effect  is  probably  not  too  significant,  since  data  by  Mentzoni  indicate  a 
maximum' factor-of -two  effect  at  relatively  low  temperature. 

3. 1.4. 3  Pulsed  Signal  Breakdown 

Flow  field  ionization  can  affect  pulsed  signal  breakdown  in  terms  of  the  separation  constant 
y  .  which  is  logarithmic  in  the  initial  electron  density.  The  relative  importance  of  this 
effect  is  somewhat  difficult  to  define,  since  it  depends  on  two  arbitrary  parameters,  the 
initial  and  final  electron  density.  The  final  electron  density  level  usually  used  in  the 
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_o 

breakdown  criterion  is  the  critical  electron  density,  which  in  units  of  cm  is  approxi- 
-8 

mately  equal  to  10  times  the  square  of  the  microwave  frequency  (for  example,  see 
References  14  and  15).  Most  authors  generally  state  that  when  the  initial  ionization  from 
the  flow  field  is  already  equal  to  or  greater  than  the  critical  value,  then  the  breakdown 
field  is  essentially  zero,  since  the  gas  may  already  be  considered  to  be  broken  down  and 
there  can  be  no  penetration  of  the  microwave  field  through  the  plasma. 

On  the  other  hand,  the  CW  breakdown  criterion  normally  used  does  not  depend  on  the 
electron  density  level,  but  only  on  the  rate  of  increase  of  electron  density  from  the  micro- 
wave  field.  Thus  no  such  critical  electron  density  effect  as  above  appears.  But  the  two 
physical  situations  cannot  be  this  different  merely  because  one  is  for  pulsed  signals  and 
the  other  is  for  CW.  The  difficulty  is  obviously  a  result  of  the  arbitrary  criteria  used  in 
defining  breakdown. 

In  the  case  of  the  reentry  induced  flow  field  of  a  slender  cone  body,  it  is  possible  to  have 

greater  than  critivid  electron  density  in  a  thin  layer,  such  that  total  signal  blackout  does 

not  occur.  Thus  there  will  be  a  significant  transmission  of  power  through  the  plasma  which 

win  increase  linearly  with  input  power  until  breakdown  is  approached,  when  the  transmitted 

power  will  attenuate  nonlinearly  due  to  the  thickening  of  the  plasma  layer  by  the  breakdown 

process.  This  indicates  that  the  only  physically  meaningful  way  to  calculate  breakdown  in 

the  presence  of  strong  flow  field  ionization  is  to  calculate  signal  transmission  as  a  function 

of  input  power  level,  taking  into  account  the  increased  level  of  ionization  induced  by  the 
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field,  as  has  been  suggested  by  Epstein  . 

A  careful  consideration  of  the  form  of  the  breakdown  criterion  is  not  only  necessary  for  pulsed 
signal  breakdown  in  the  presence  of  ionization  but  also  for  breakdown  with  convection,  as 
pointed  out  in  the  discussion  at  the  end  of  Section  3. 1. 2.  A  meaningful  prediction  of  break¬ 
down  in  the  presence  of  flow  field  ionization  with  convection  requires  one  to  calculate  the 
electron  density  distribution  over  the  antenna  and  the  power  transmitted  through  this  plasma 
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as  functions  of  the  input  power.  The  importance  of  this  consideration  lies  mainly  in  its 
effect  on  the  breakdown  criterion  and  not  so  much  in  the  phenomenon  of  alteration  of  the 
antenna  near-field  and  impedance  by  the  plasma,  although  the  latter  may  be  important 
considerations  if  the  electron  density  is  initially  very  high. 

3.2  THE  NET  IONIZATION  FREQUENCY 

The  net  ionization  frequency  is  the  most  important  term  in  the  continuity  equation  for 
microwave  breakdown.  This  quantity,  through  its  dependence  on  the  microwave  field 
intensity,  determines  the  breakdown  field.  In  practice  the  breakdown  field  is  obtained  by 
first  solving  the  continuity  equation  for  the  net  ionization  frequene*  under  the  appropriate 
breakdown  criterion.  Then  the  dependence  of  the  net  ionization  frequency  on  the  electric 
field  determines  the  breakdown  field. 


3.2.1  DEPENDENCE  ON  MICROWAVE  FIELD 

Unfortunately,  the  ionization  frequency  has  never  been  measured  directly  except  for  hydrogen, 
for  which  it  was  done  by  an  optical  method  by  Cottingham  and  Buchsbaum  .  The  data  for 
air  have  been  inferred  from  either  diffusion  controlled  or  pulse  controlled  breakdown  fields, 
using  some  form  of  phenomenological  theory. 
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Two  representations  of  the  data  for  air  exist.  MacDonald  has  reduced  the  data  of  Herlin 
24 

and  Brown  to  the  form  of  curves  of  vRA  versus  pA.  with  EA.  as  a  parameter,  where  A.  is  the 

wavelength  and  E  is  the  electric  field.  In  the  MacDonald  method  of  predicting  breakdown, 

one  uses  these  curves  in  conjunction  with  Equation  3.1.4.  2-10,  using  the  following  expression 

for  the  effective  field  E  : 

e 


25 

Scharfman  and  Morita  have  obtained  data  by  a  focused  pulsed  microwave  technique.  They 
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have  interpreted  their  data  in  terms  of  a  single  curve  of  net  ionization  frequency  as  a 
function  of  E^/p,  including  in  the  results  the  data  of  several  other  investigators  also  reduced 
to  this  form.  Their  results  are  shown  in  Figure  3. 2. 1-1  along  with  a  simple  power  law 
approximation  given  by  the  following  equation: 


(3.  2. 1-2) 


This  function  has  advantages  in  terms  of  nonuniform  electric  field  problems  (see  below). 

3.  2.  2  HIGH  GAS  TEMPERATURE 

At  the  present  time  the  effect  of  high  gas  temperature  on  the  net  ionization  frequency  is 
not  well  understood  theoretically,  but  there  is  enough  experimental  data  to  indicate  the 
existence  of  an  effect  in  shock  heated  air  (see  References  16,  26,  27,  and  28).  Some  of 
the  more  pertinei  data  are  shown  in  Figure  3.  2.  2-1. 

16 

Light  and  Taylor  have  proposed  a  model  for  the  high  temperature  effect  on  the  ionization 
frequency  which  relates  phenomenologically  to  the  dependence  of  the  electros  energy  distri¬ 
bution  on  vibrational  excitation  of  oxygen  molecules.  Their  equation  is  as  follows: 


(3. 2. 2-1) 


where  T  is  the  gas  temperature  in  °K,  X/Xq  is  the  ratio  of  the  mole  fraction  of 
oxygen  molecules  to  their  mole  fraction  in  cold  air  and  F  [(Ee/p*)2]is  the  net  ionization 
frequency  function  in  cold  air.  The  above  equation  represents  the  apparent  trends  of  the 
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data  fairly  well  if  the  dissociation  of  oxygen  is  relatively  complete*  at  the  higher  tempera¬ 
tures  in  Light  and  Taylor's  experiment;  however,  the  range  of  temperatures  and  field 
strengths  is  rather  small  for  the  high  temperature  data. 

3.  2. 3  CHEMICAL  COMPOSITION  EFFECTS 

It  is  obvious  that  the  net  ionization  frequency  as  a  function  of  electric  field  (and  temperature) 
must  depend  on  the  gas  composition,  for  several  reasons.  The  ionization  potentials  of 
many  gases  are  lower  than  air  and  there  are  also  possibilities  of  many  additional  mechanisms 
for  energy  and  charge  transfer,  all  of  which  affect  the  ionization  frequency.  The  attachment 
frequency  can  be  affected  very  significantly  also  with  either  the  addition  of  electrophilic 
materials  acting  to  enhance  it  or  the  dissociation  of  oxygen  molecules  acting  as  a  damper. 

g 

(See  the  previous  report  in  this  series  for  a  brief  discus rion  of  some  of  the  more  recent 
data  on  these  effects. ) 

3.3  GRADIENTS  OF  ANTENNA  FIELDS  AND  FLOW  FIELD  PROPERTIES 
The  continuity  equation  discussed  in  detail  above  is,  in  general,  an  elliptic  partial  differential 
equation  of  second  order  involving  both  time  and  space  derivatives  and  possessing  nonconstant 
coefficients.  It  is  apparent  from  this  fact  and  from  some  of  the  discussions  above  that  the 
mathematical  criterion  used  in  defining  breakdown  in  terms  of  a  solution  of  this  equation  will 
have  a  significant  influence  on  the  results.  Moreover,  the  physical  effects  of  gradients  of 
the  antenna  near-field  and  the  flow-field  properties  will  be  significant  in  the  problem  of 
reentry  antenna  breakdown.  Thus  these  gradients  must  be  taken  into  account  in  solving 
the  continuity  equation  to  predict  breakdown. 

3.3.1  ANTENNA  NEAR -FIELDS 

The  problem  of  nonuniform  antenna  near-fields  has  been  discussed  by  Brown  and  MacDonald29 
in  terms  of  a  limiting  boundary  on  the  region  of  validity  of  the  diffusion  theory  of  breakdown. 

*  The  thermodynamic  equilibrium  concentrations  of  02  are  needed  to  explain  the  data. 
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Within  this  region  of  validity  the  continuity  equation  (in  the  absence  of  convection)  supposedly 

can  be  solved  as  a  homogeneous  differential  equation  with  constant  coefficient,  giving  as  a 

2 

solution  the  eigenvalue  (  -1/A  ),  where  A  is  the  diffusion  length.  Brown  and  MacDonald 
assume  that  the  uniform  field  assumption  is  correct  in  the  parallel  plate  situation  if  the 
plate  separation  is  less  than  or  equal  to  a  half  wavelength,  so  that  no  more  than  a  single 
loop  of  a  standing  wave  of  electric  field  would  be  maintained.  This  gives  the  result  that  the 
uniform  field  assumption  works  if  A  is  less  than  X  /  (2  7r). 

30  31  32 

Herlin  and  Brown  *  and  MacDonald  have  calculated  breakdown  in  nonuniform  fields  by 

33 

a  method  which  is  discussed  fully  by  MacDonald  ,  using  the  following  form  of  continuity 
equation: 

+(~d)q  (^")/3  *  =  °  (3.3.1-1) 

where  0  is  held  constant  throughout  the  gas  volume.  The  results  for  cylindrical  cavities 
depend  on  the  geometry  and  they  indicate  that  the  effective  height  of  the  cavity  can  be  much 
less  than  its  actual  height  if  the  radius  is  relatively  small.  For  a  spherical  cavity  the 
effective  size  of  the  cavity  can  be  reduced  by  a  factor  of  two  or  three  by  this  effect.  Of 
course  for  the  case  of  antenna  breakdown,  boundaries  such  as  exist  in  cavities  are  absent 
and  one  must  solve  an  equation  similar  to  Equation  3.3. 1-1  for  the  conditions  imposed  by 
the  antenna  near-field.  Most  investigators  of  antenna  breakdown  assume  a  diffusion  length 
either  by  analogy  with  parallel  plate  breakdown  or  by  fitting  experimental  breakdown  data  to 
diffusion  theory. 

The  determination  of  the  magnitude  of  the  antenna  near-fields  is  necessary  for  prediction 
of  the  breakdown  condition  shown  in  Equation  3. 3. 1-1.  We  consider  three  antennas: 


a.  Electric  dipole 

b.  Magnetic  dipole  (electric  current  loop) 

c.  Conical  equiangular-spiral  antenna  (C ESA) 


The  electric  dipole  has  been  used  in  previous  analyses34  and  is  a  useful  approximation  to 

several  radiating  configurations.  The  magnetic  dipole  is  used  in  the  analyses  of  slots  on 
35 

planes  .  The  CESA  is  of  interest  in  ECM  applications  for  its  broadband  characteristics 
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and  practical  geometry  ,  although  polarization  diversity  may  be  a  problem3^. 

The  analysis  of  the  fields  of  electric  and  magnetic  dipoles  provides  useful  and  simple 

approximations  to  the  near-fields  of  more  complicated  antennas.  As  indicated  in  the  first 
38  — * 

quarterly  report  ,  the  electric  field  E,  is  found  from  the  vector  potential,  A: 


-J-  v  .  v 


(3.  3. 1-2) 


with 


J  (r',  O',  0') 


dr* 


where 


(3.  3. 1-3) 


j  (I-*,  e»,  0*) 

-  Antenna  source  current  density 

K 

=  27rA 

\ 

=  Radiated  wavelength 

r*,  0',  0* 

—  Spherical  components  of  source  current  elements 

dr* 

=  Source  volume  element 

— *  — ♦  I 

r  -  r*  |  =  r 

f  2  2  2  1 

=  [<X  -  X')  +  (Y  -  Y »)  +  (Z  -  Z»)2  J 

For  the  electric  dipole  (Figure  3. 3. 1-1)  we  use  as  a  source  current  I  and  we  have 
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(3.  3. 1-5) 


£  t  i.  „  “J  Kr 


A  =  KA  =  - ---  6 
z  4nr 


where  &  is  the  length  of  the  dipole. 


-  A  A 

E  =  r  E  +0  Efl 
r  0 


(3.  3. 1-6) 


r  _  M  -JKr  /  Z  1  \ 

Er'47re  (~2  +  r)COS0 

\r  jwer  / 

\  r  jwer  / 


(3.  3.1-7) 


(3.  3. 1-8) 


where  Z  =^z/e.  Plots  of  Ef,  EQ,  and  E  »  Je  2  +  " 
shown  in  Figures  3.  3. 1-2  through  3.  3. 1-7. 


q  versus  r  and  0  in  the  near-field  are 


The  magnetic  dipole  (see  Figure  3.  3. 1-8)  is  characterized  by  its  strength  KS  =  I  ir  a2. 


A  =  0  A 


(3.  3- 1-9) 


A  _  Iffa 
0  47T 


a_  /IK  ,  1  \  -jKr  . 

•  +  ~2 ) e  sin6» 


(3.3.1-10) 


E  =  0  E 


=  ^a2  /K2  _  ]K_\  -jKr 


0  4tt  ^  r  ^2 


sin  0 


(3.3.1-11) 
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Figure  3.  3. 1-2.  Magnitude  of  Electric  Field  of  Short  Electric  Dipole,  Normalized  to  Magnitude  of  Field  at 

R  =  x/T  |E|  =  / E  2  +  ;  Axial  Radiation 


Expanded  View  of  F 


Figure  3.  3. 1-4.  Magnitude  of  Electric  Field  of  Short  Electric  Dipole,  Normalized  to  Magnitude  of  Field  at 

R  =  \/2~,  j  E  |  =  / E  2  +  Efl2;  Broadside  Radiation 
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Figure  3.  3. 1-6.  Polar  Plot  of  |E  |  for  Short  Electric  Dipole,  Normalized  to  |E  |  at  R  = 
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Figure  3. 3. 1-8.  Configuration  of  Magnetic  Dipole  (Equivalent  to  a  Small  Electric  Current 
Loop)  Relative  to  the  Electric  Dipole  of  Figure  3. 3. 1-1 
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The  field  strength  vs  r  and  0  is  plotted  in  Figures  ”.  3. 1-9  through  3. 3. 1-17. 

The  CESA  has  been  analyzed  theoretically  by  solving  the  integral  Equation  3.  3. 1-3  for  the 
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antenna  current  distribution  .  The  radiation  patterns  are  well  known.  Extensive  measure¬ 
ments  have  been  made^’  ^  of  the  radiation  patterns  and  near-fields.  These  are  reproduced 
in  Figures  3.  3. 1-18  and  3.  3. 1-19.  The  antenna  current  J  in  general  has  three  components, 
so  that  the  near-fields 

J  =  rJ  +  0  J.  +  0  J_  (3»  1_12) 

r  0  0 

will  contain  both  radial  and  azimuthal  nomniformities.  The  field  expressions  are  thus 

40 

complicated  infinite  series  in  terms  of  Bessel  and  Legendre  functions  . 

A  practical  approach  to  obtaining  approximate  values  for  the  diffusion  length  of  typical 

41 

antennas  might  utilize  measured  values  of  the  near-field  in  the  following  way.  For  a  10 

degree  conical  spiral,  the  width  of  the  active  region  is  approximately  0. 45  X  if  the  length  of 

the  near-field  is  taken  between  10  db  points  of  Figure  3.  3. 1-18.  The  impedance  measure- 

36 

ments  of  Figure  3.  3. 1-19  indicate  a  10  to  1  bandwidth  is  possible.  Other  investigators 
have  not  achieved  such  wide  bandwidth,  possibly  due  to  their  narrow  band  balun.  The  active 
region  shown  in  Figure  3.  3. 1-18  is  the  effective  radiating  aperture  of  the  antenna.  That  is, 
although  the  CESA  is  an  end  fire  antenna  as  shown  in  Figure  3.  3. 1-19,  the  peak  fields  occur 
in  the  active  region.  We  therefore  make  the  approximation  that  the  CESA  may  be  treated  as 
a  slot  of  width  b  where  we  equate  b  to  the  length  of  the  active  region.  Using  this  approxima¬ 
tion  we  may  obtain  an  order  of  magnitude  diffusion  length  for  the  conical  equiangular  spiral, 
by  using  the  well  known  result  for  a  slot  of  width  b,  namely  A=  b/(2ir).  Thus  we  estimate 
X  =  14A. 
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Magnitude  of  Electric  Field  of  a  Small  Current  Loop,  Normalized  to  Magnitude  of  Field  at 

R  =X/50;  Oblique  Radiation 
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Figure  3.  3. 1-18.  Polar  Plots  of  |E^|for  Small  Electric  Current  Loop,  R  =  \/2,  3X/5, 


3.  3.  2  FLOW  FIELD  GRADIENTS 

The  above  discussion  indicates  the  importance  of  field  gradients  in  antenna  breakdown.  It 
is  possible  to  consider  such  gradients  as  equivalent  to  an  additional  electron  loss  mechanism. 
In  a  similar  way,  the  nonuniform  flow  field  around  a  reentry  vehicle  can  be  effectively  an 
electron  loss  mechanism.  Among  these  effects  are  gradients  of  gas  density,  convective 
velocity,  temperature,  composition,  and  electron  density.  Thus,  all  the  terms  in  the 
coefficient  of  iii  in  the  continuity  equation  are  subject  to  gradients.  The  complete  solution  of 
the  reentry  breakdown  problem  must  consider  these  nonuniform  properties,  as  Fante*°  has 
indicated  in  an  approximate  analysis  using  a  two-layered  model. 

3.4  ILLUSTRATIVE  CALCULATIONS 

The  continuity  equation  ir  t.he  reentry  environment  can  be  written  in  the  following  form  as  a 
breakdown  condition*: 


e 


assuming  an  effective  diffusion  length  A  ,  transitional  diffusion  coefficient  D  ,  effective 

e  s 

velocity  V,  and  initial  flow  fielc.  electron  density  n^  can  be  defined.  A  convenient  set  of 
proper  variables  is  (p*X),  (X/A),  and  the  number  N  of  rf  cycles  in  a  pulse  of  duration  t. 


N  (p*\) 


*<?) 


V 


4  (DgP*) 


(3.4-2) 


Now  the  reentry  effects  can  be  written  in  terms  of  a  number  of  parameters,  including  (p/p*), 

(D  /D  ),  (A/A  ),  y  ,  andy, ,  where  A  is  the  diffusion  length  of  the  antenna  in  the  absence  of 
s  -  e  o  1 

flow  field  gradients. 

*As  pointed  out  above,  this  form  of  velocity  effect  may  be. conservatively  large,  but  it  is 
easily  changed  by  dropping  the  last  term  and  changing  t  to  L/V,  whichever  is  smaller. 
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3. 1-20.  Amplitude  of  Near-Fields  Measured  With  a  Small  Shielded  Loop  Probe  Along  the  Surface  of  One 
Electric  Far-Field  Radiation  Patterns  Corresponding  to  a  Truncation  at  Indicated  Points  (From  J.  D. 
Dyson,  Trans.  IEEE,  Ant.  and  Prop.  ,  Vol  AP-15,  Page  488,  July  1965) 


c 


(3.4-3) 


(3.4-4) 


where  c  is  the  speed  of  light,  p  is  the  free  stream  pressure  (p  -  p*),  and  nQ  is  the  initial 
electron  density  in  the  ordinary  pulse  breakdown  criterion. 


(3. 4-5) 


A  littlo  manipulation  gives  the  result  in  the  following  form: 


(3,4-6) 


where  the  CW  breakdown  ionization  frequency  without  reentry  effects  is 


(3.4-7) 


This  equation  shows  the  effects  of  reentry  directly  in  each  of  its  terms.  The  first  term 
shows  that  diffusion  is  affected  by  the  transitional  diffusion  coefficient  and  by  gradients. 

The  second  term  shows  that  the  convection  effect  is  important  at  high  pressure  and  when 
diffusion  is  inefficient.  The  last  two  terms  show  that  flow  field  ionization  tends  to  decrease 
the  effect  of  pulsed  signals. 
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The  effects  of  reentry  on  the  breakdown  power  of  an  antenna  can  be  related  to  the  ionization 

frequency  by  using  Equation  3. 2. 1-2.  Consider  a  simple  illustrative  case,  using \/A=  2tt, 

3  8  d 

N  =  10  ,  n  /n  =  n  /n  10  ,  p*/p  =  4.  5,  v  =  2  :c  10  feet/second,  and  A/A  =3.  In  the 
•D  o  i  o  e 

presence  of  ionization  we  reduce  to  0.  01  D  .  (Note  that  in  the  finite  time  theory  of 
convection  the  pulse  length  used  in  this  example  is  equivalent  to  an  antenna  length  of  0.  02  ) 

Table  3.4-1  shows  the  results  for  p  X  =  36  cm  torr. 

These  results  show  that  the  effect  of  convection  is  weak  in  the  absence  of  ionization,  except 
in  the  space  charge  coupled  theory,  in  which  case  it  represents  the  largest  effect  shown. 

The  latter  theory  is  not  shown  for  the  case  of  ionization  (ambipolar  diffusion),  since  it  is  then 
identical  with  the  conventional  theory. 

In  addition  to  the  effects  shown  in  the  table,  the  effects  of  gas  temperature  and  composition 
on  the  net  ionization  frequency  must  be  considered.  As  shown  in  Figure  3. 2. 2-1,  the  gas 
temperature  effect  can  be  on  the  order  of  6  db  in  the  shock  layer.  At  any  rate,  the  gas 
temperature  effect  will  iend  to  reduce  the  numbers  shown  in  the  last  column  of  Table  3. 4-1. 
Gas  composition  may  have  either  kind  of  effect. 

This  type  of  table  is  useful  in  predicting  where  breakdown  will  occur,  by  searching  for  the 
lowest  db  values.  For  example,  if  the  conventional  convection  theory  is  correct,  then  CW 
breakdown  will  occur  in  the  free  stream,  assuming  that  the  antenna  field  does  not  fall  off  too 
rapidly.  Pulse  breakdown  may  be  affected  by  gas  temperature  enough  to  cause  breakdown 
in  the  shock  layer.  If  the  space  charge  coupled  theory  is  true,  then  the  presence  of  shock 
layer  ionization  (there,  conventional  convection  theory  is  equivalent  to  space  charge  coupled 
theory)  will  make  the  free  stream  and  shock  layer  roughly  equivalent.  High  gas  temperature 
will  then  cause  breakdown  in  the  shock  layer. 

The  above  discussion  applies  only  to  one  value  of  p\,  which  has  been  chosen  near  the  normal 
breakdown  minimum.  The  results  are  more  generally  useful  if  we  consider  the  net  ioniza¬ 
tion  frequency  as  a  measure  of  the  importance  of  a  given  effect,  keeping  in  mind  that,  relative 
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Table  3.4-1.  Typical  Effects  of  Reentry  on  Antenna  Breakdown 


Conditions 

V  /p*  (sec  torr)  1 
n 

Breakdown  Power,  db 

Free  stream: 

CW 

4.  9x1 04 

0 

Pulse 

1.  5xl07 

+8.1 

CW,  v 

l.lxlO5 

+0.8 

CW,  v* 

5.  9x10® 

+14.1 

Pulse,  v 

1.  5x1 07 

+8.1 

Pulse,  v* 

5.  9x1 08 

+14.1 

Shock  layer: 

CW 

1.  lxiO4 

+9.5 

CW,  n1 

1.  lxlO2 

+9.3 

CW,  n  ,  A 
l  e 

lxlO3 

+9.  3 

CW,  n  ,  A  ,  v 
l  e 

1.  3x10 

+14.2 

CW,  A  ,  v 
e 

1.  lxlO5 

+10.  8 

CW,  A  ,  V* 
e 

1. 2x10® 

r* 

+21.  5 

Pulse 

3.  4x10° 

+15.7 

Pulse,  n^ 

1.  lxlO2 

+9.3 

Pulse,  n, ,  A 

1  e 

lxlO3 

+9.3 

Pulse,  n, ,  A  ,  v 

1  e 

1. lxlO5 

+14.2 

Pulse,  A  ,  v 
e 

3.  5x1 06 

+15.7 

Pulse,  A  ,  v* 
e 

1.2x10® 

+21.5 

NOTE:  In  the  Conditions  column,  ^  refers  to  ionization  effects,  refers  to  gradients,  v 
refers  to  convection,  and  v*  is  convection  in  the  space  charge  coupled  theory. 


2 

to  the  diffusion  effect,  the  convection  effect  scales  as  (p\)  and  the  pulse  effect  as  p\.  For 
example,  on  the  low  pressure  side  of  the  breakdown  curve  the  breakdown  power  is  proportional 
to  the  square  of  the  net  ionization  frequency  and  there  is  a  13  db  increase  in  free  stream 
breakdown  field  over  the  shock  layer  values.  It  is  easily  seen  then  that  at  one-tenth  of  the 
pressure  discussed  above,  the  diffusion  effect  would  dominate  over  all  others  and  the  only 
significant  reentry  effects  would  be  ionization  and  gradients. 

3.5  RECOMMENDATIONS 

From  the  above  discussions  it  is  evident  that  the  theoretical  prediction  of  antenna  breakdown 
in  a  reentry  environment  presents  a  number  of  problems  at  the  present  time.  The  most 
serious  theoretical  problem  is  the  effect  of  high  gas  temperature,  where  theoretical  models 
are  practically  nonexistent  and  experimental  data  cover  only  a  narrow  range  of  conditions. 
Therefore  effort  is  needed  in  the  following  areas  of  development: 

a.  Breakdown  data  at  high  temperature  is  needed  over  a  wider  range  of  laboratory 
experiment  conditions. 

b.  Phenomenological  models  need  to  be  developed  and  checked  against  the  data. 

c.  Crucial  experiments  are  needed  to  verify  the  models  once  they  have  been  derived. 

The  convective  effect  is  important  only  in  view  of  the  widely  different  theoretical  approaches 
presently  extant.  The  older  theories  can  be  shown  to  be  approximately  equivalent  in  the 
broad  sense,  and  on  the  basis  of  these  alone  the  effect  of  convection  would  not  be  significant 
in  reentry.  However,  the  coupled  convection-diffusion-space  charge  theory  predicts  much 
stronger  effects.  More  theoretical  work  is  needed  in  formulating  breakdown  criteria  and 
resolving  inconsistencies  among  the  theories.  Laboratory  experiments  over  an  extended 
range  of  conditions  of  velocity,  electron  density,  and  diffusion  length  are  needed  to  settle  the 
observational  facts  and  to  give  crucial  tests  of  the  theory.  At  the  present  time  it  appears 
that  some  experimenters  see  one  effect  and  others  see  the  other. 
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Careful  theoretical  analyses  are  needed  in  order  to  predict  the  breakdown  of  an  antenna  in 
the  reentry  environment  in  terms  of  a  solution  of  the  continuity  equation.  The  factors  which 
need  careful  attention  are  breakdown  criteria  in  terms  of  effects  of  ionization  and  convection, 
near-field  gradients,  and  gradients  of  flow  field  properties. 

The  fact  that  many  flow  field  parameters  in  addition  to  ionization  tend  to  affect  breakdown 
means  that  the  effects  of  reentry  will  be  difficult  to  isolate  in  a  flight  experiment.  In  view 
of  the  fact  that  flight  experiments  should  be  conducted  as  a  final  step  in  verifying  the 
importance  of  various  reentry  effects,  it  is  a  necessity  to  take  all  possible  steps  to  design 
such  flight  experiments  to  give  the  needed  isolation  of  the  various  important  effects  of 
reentry.  This  consideration  is  just  as  important  as  doing  a  straightforward  system  test 
experiment  by  flying  a  given  antenna  and  vehicle  system,  since  the  system  test  results  are 
difficult  to  extrapolate  to  other  situations. 
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SECTION  4 

EXPERIMENT  DESIGN  REQUIREMENTS 

This  section  is  intended  to  give  the  general  relationships  between  the  design  requirements 
and  the  objectives  of  antenna  breakdown  experiments.  It  derives  requirements  on  the  signal 
characteristics  (pulse  length)  and  the  field  geometry  (antenna  design)  needed  to  meet  various 
experiment  objectives.  A  decision  to  design  an  experiment  for  a  given  set  of  detailed 
objectives  would  then  lead  to  specific  requirements  on  the  pulse  length  and  antenna  design. 
Requirements  on  pulse  amplitude  automatically  result  from  these  two  requirements. 

Various  possible  objectives  of  antenna  breakdown  experiments  can  be  anticipated,  but  these 
are  generally  of  two  types,  namely  system  tests  and  basic  data  experiments.  The  objective 
of  a  system  test  is  to  obtain  data  on  the  extent  to  which  the  reentry  environment  changes  the 
antenna  power  handling  capability  relative  to  that  in  a  static  environment,  for  a  given  antenna, 
signal,  and  vehicle  system.  This  type  of  experiment  gives  only  a  single  "system  data  point" 
in  that  its  results  cannot  rightfully  be  extrapolated  to  other  systems.  On  the  other  hand,  the 
general  objective  of  a  basic  data  type  of  experiment  is  to  yield  data  which  can  be  extrapolated 
to  other  systems,  because  it  gives  basic  data  on  the  magnitude  and  form  of  the  effect  which 
a  given  specific  feature  of  the  reentry  environment  has  on  antenna  breakdown. 

The  technical  requirement  of  a  system  test  type  of  experiment  is  that  it  simulate  all  of  the 
reentry  effects  without  compromising  the  ability  to  obtain  usable  data.  The  technical 
requirement  of  a  basic  data  experiment  is  that  it  provide  data  on  the  effect  of  only  one  of  the 
reentry  environmental  factors  at  a  time. 

The  basis  for  the  establishment  of  experiment  requirements  can  be  derived  from  the  equations 
of  Section  3. 4,  in  which  the  quantities  N  and  \/A  enter  as  controllable  parameters.  N  is 
proportional  to  the  pulse  length  and  within  limits  \/A  is  ?>  function  of  the  antenna  geometry. 
The  object  of  this  section  is  to  obtain  the  condition  on  these  two  parameters  for  which  each 
of  the  various  reentry  effects  dominates.  The  result  will  be  expressible  in  terms  of  a  graph 
of  N  versus  \/A  showing  regions  in  which  each  of  these  effects  dominates. 
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The  equation  which  we  use  as  the  basis  for  these  calculations  is 


(4-1) 


where  v,  and  v  are  the  ionization  and  attachment  frequencies,  respectively.  In  order  to 

1  a  23  4 

give  realistic  representation  of  the  attachment  we  use  MacDonald's  value  of  2. 12  x  10 

(sec  torr)  1  for  v  /p. 

Si 


Before  proceeding  to  the  specification  of  requirements  on  the  controllable  parameters  N  and 
v/A,  it  is  helpful  to  take  account  of  the  fact  that  transmitter  power  output  is  always  a 
practical  consideration.  Thus  it  is  necessary  to  conduct  experiments  in  an  altitude  regime 
near  the  breakdown  minimum,  which  we  assume  to  correspond  to  a  value  of  36  cm  torr  for 
pX.  The  boundaries  of  the  regions  in  which  each  of  the  terms  in  this  equation  dominates  over 
the  sum  of  all  of  the  others  are  then  given  by  simple  equations  between  N  and  X/A,  which 
are  plotted  in  Figure  4-1.  In  this  figure  the  regions  of  the  N  versus  X/A  space  marked  by 
arrows  are  dominated  by  the  different  effects  indicated.  In  general,  the  shapes  of  these 
boundaries  will  be  the  same  regardless  of  what  part  of  the  breakdown  curve  is  considered. 

As  is  evident,  the  value  of  the  diffusion  coefficient  has  a  strong  effect.  Note  also  that  the 
slight  separation  between  the  convection  boundary  and  the  pulse  and  diffusion  boundaries  is 
caused  by  attachment  being  nearly  equal  to  convection  in  the  free  diffusion  regime. 


It  is  possible  to  construct  similar  boundaries  for  th  ;  breakdown  control  regions  for  other 
values  of  velocity  arid  pX,  but  it  is  more  convenient  to  consider  only  the  asymptotes  of  these 
curves,  keeping  in  mind  the  general  shapes  as  modeled  by  Figure  4-1.  Different  conditions 
merely  move  these  asymptotes  around  the  graph  without  rotation. 
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The  pulse-diffusion  control  boundary  has  the  asymptotic  form  of 


(4-2) 


The  convection  effect  dominates  over  attachment  at  velocities  of  12,000  and  1200  feet/ 
second  for  free  and  ambipolar  diffusion,  respectively.  (The  space  charge  coupled  theory 
would  give  a  value  of  120  feet/second  for  free  diffusion. )  The  diffusion-convection  boundary 
is 


\  =  V  (M) 
A  2  (D  p) 


=  3. 1  x  10"7  V  (p\) 


and  the  diffusion- attachment  boundary  is 


The  pulse-convection  boundary  is 


(4-3) 


(4-4) 
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N 


4c  D  p 
(pX)  v2 


3.  5  x  10 

(pX)  v2 


and  the  pulse-attachment  boundary  is 


WV/p)  fa(nQ  ) 


2.<;x  io7 

(pX) 


(4-5) 


(4-6) 


Note  that  the  space  charge  coupled  convection  theory  reduces  N  at  the  pulse-convection 
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boundary  by  a  factor  of  10  . 


In  order  to  extend  the  conditions  for  these  criteria  beyond  the  breakdown  minimum,  consider 
the  high  pressure  side  of  the  curve,  which  is  dominated  by  the  effective  field  effect: 


(pX)  -  36^2p/p°  -  1  {4-7) 

where  p/pQ  is  the  breakdown  power  relative  to  the  minimum  value.  Thus,  for  example, 

the  value  of  pX  is  increased  by  4.  35  when  the  breakdown  power  is  10  db  above  the  minimum 

2 

on  the  high  pressure  side,  increasing  N  (X/A)  at  the  pulse-diffusion  boundary  by  the  same 
factor.  Also  X/A  at  the  attachment  or  convection-diffusion  boundary  is  increased  by  the 
same  amount.  The  value  of  N  at  the  attachment  or  convection-pulse  boundary  is  reduced 
by  this  factor. 


4-5 


On  the  low  pressure  side  of  the  curve,  only  diffusion  or  pulse  effects  can  be  important  and 
the  attachment/convection  region  does  not  exist.  The  effective  field  on  this  side  of  the 
curve  is  given  by 


E 

e  EX 

p  36 


while  at  the  minimum  breakdown  field, 


Thus,  if  we  use  the  power  law  for  the  ionization  frequency, 


Simultaneous  solution  of  Equations  4-2  and  4-10  gives 


(4-8) 


(4-9) 


(4-10) 


(4-11) 


The  pulse  diffusion  boundary  at  10  db  up  from  the  minimum  is  then  given  by 


(4-12) 


If  we  consider  5  and  2  as  the  smallest  practical  values  of  N  and  A/A,  respectively,  then 
solving  for  E/Eq,  we  find  that  for  free  diffusion  pulse,  controlled  breakdown  can  be 
maintained  only  up  to  18  db  above  the  minimum  (8  times  the  minimum  field)  on  the  low 
pressure  side.  Similar  practical  limitations  appear  in  other  factors  involving  diffusion 
controlled  or  pulse  controlled  breakdown  on  the  high  pressure  side.  However,  it  is 

gt  ne rally  true  that  an  approach  like  that  given  above  is  necessary  and  useful  to  the  design 
of  meaningful  experiments. 


The  following  general  considerations  enter  into  the  generation  of  experimental  requirements 
on  pulse  length  and  antenna  design.  As  indicated  in  Figure  4-2,  an  experiment  designed  to 
isolate  one  of  the  six  phenomena  listed  requires  the  experiment  design  to  lie  in  the  region 
indicated.  Regions  2  and  3  are  pulse  and  diffusion  regions,  respectively,  and  region  1  is 
either  convection  or  attachment.  (Gas  density  experiments  are  less  interesting  than  the 
others,  because  this  effect  is  well  understood. )  A  gas  temperature  experiment  would  best 
be  done  in  region  2.  Ionization  experiments  could  be  done  in  either  region  3,  using 
diffusion  as  an  indicator,  or  region  1,  using  convection  as  an  indicator.  Effects  of  additives 
on  attachment  and  ionization  frequency  could  be  isolated  by  doing  an  experiment  with  both 
long  and  short  pulses,  provided  convection  were  eliminated  from  region  1.  The  short  pulse 
would  isolate  the  effects  on  the  ionization  frequency  and  the  long  pulse  would  give  the  addi¬ 
tional  effect  ce  attachment.  Gas  temperature  would  have  to  be  either  minimized  or  accounted 
for  by  "turning  off’  the  additive  injection. 

The  design  of  a  system  test  experiment  is  made  more  meaningful  by  the  above  techniques 
also.  Generally  it  is  necessary  to  conduct  antenna  breakdown  experiments  with  a  train  of 
pulses  of  various  amplitudes,  the  pulses  which  break  down  indicating  the  breakdown  power 
level.  But  care  must  be  taken  so  as  not  to  change  the  effects  of  reentry  by  using  pulses 
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EXPERIMENT 


REQUIRED  REGION 


CONVECTION 

DENSITY 

GRADIENTS 

TEMPERATURE 

IONIZATION 

ADDITIVES 

1  =  CONVECTION  REGION 

1*  =  ATTACHMENT  REGION 


1 

2 

3 

2 

3  (OR  1) 
1*  (AND  2) 


Figure  4-2.  Pulse  Length  and  Antenna  Geometry  Requirements  for  Basic  Data  Experiments 
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which  are  too  short.  Thus  the  experimental  system  must  lie  on  the  same  side  of  the  pulse 
control  boundary  as  the  actual  system  being  tested.  If  the  system  lies  in  region  2,  then 
the  experiment  must  be  designed  with  the  same  pulse  length  as  the  system  itself. 

To  summarize  the  results  of  this  section,  we  tabulate  numerical  values  for  the  boundaries 
of  the  various  control  regions,  using  a  flow  velocity  of  20, 000  feet/second.  For  the 
diffusion-convection/attachment  boundary  we  get  the  following  results  for  the  value  of  A /A 
above  which  diffusion  dominates: 


Breakdown 

Minimum 

10  db  Above  Minimum 
at  High  Pressure 

Convection,  with  free  diffusion 

6.8 

30 

Convection,  with  ambipolar  diffusion, 
or  with  free  diffusion  in  the  space 
charge  coupled  theory 

680 

3000 

Attachment,  with  free  diffusion 

4.1 

18 

Attachment,  with  ambipolar  diffusion 

41 

180 

For  the  pulse-convection/attachment  boundary, 

dominates  is  as  follows: 

the  value  of  N  below  which  the  pulse  length 

Breakdown  10  db  Above  Minimum 
Minimum  at  High  Pressure 


Convection,  with  free  diffusion 

2.6xl04 

_ ~  * - -  w 

6xl03 

Convection,  with  ambipolar  diffusion 

260 

60 

Convection,  with  free  diffusion  in  the 
space  charge  coupled  theory 

2.6 

0.6 

Attachment 

7.2xl05 

1.7xl05 
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For  the  diffusion -pulse  boundary  the  value  of  N  (X/A)  above  which  diffusion  dominates  is 
as  follows: 

7 

Breakdown  minimum:  1.2x10 

7 

10  db  above  minimum  at  high  pressure:  5.4x10 

3 

10  db  above  minimum  at  low  pressure:  8.  5x10 

The  values  of  these  last  three  results  are  increased  by  a  factor  of  100  by  ambipolar  diffusion. 

The  following  list  summarizes  the  values  of  N  and  X/A  needed  to  maintain  the  various  control 
regions  discussed  above. 


N 

X/A 

Convection 

>  2. 6x1 04 

<6.8 

Pulse 

<  60* 

<92^ 

Diffusion** 

>  6xl05 

>  30  1 

♦Much  lower  values  are  required  in  the  space  charge  coupled  convection  theory. 

♦♦This  applies  only  to  free  diffusion.  Much  higner  values  apply  with  ambipolar  diffusion  or 
in  the  space  charge  coupled  convection  theory. 
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SECTION  5 
CONCLUSIONS 

The  vibrational  nonequilibrium  effect  on  the  flow  field  of  a  hypersonic,  slightly  blunt  body  has 
been  found  to  be  significant  at  250, 000  feet  altitude.  This  conclusion  has  been  reached  through 
results  of  calculations  using  a  superpositional  technique.  It  is  found  that  once  the  stagnation 
region  streamlines  expand  around  the  conical  region  of  the  body,  the  vibrational  temperature 
tends  to  freeze  and  approaches  the  equilibrium  state  only  after  a  long  distance.  The  vibration¬ 
al  temperature  tends  to  be  higher  than  the  equilibrium  temperature  near  the  body  and  vice 

versa,  away  from  the  body.  If  high  vibrational  temperature  decreases  the  breakdown  field 

16 

as  suggested  by  Light  and  Taylor  ,  then  vibrational  nonequilibrium  may  have  an  important 
effect  on  antenna  breakdown,  especially  near  the  vehicle  surface.  However,  boundary  layer 
effects  may  change  this  conclusion. 

Sufficient  knowledge  presently  exists  on  the  phenomena  of  reentry  which  affect  antenna  break¬ 
down  to  permit  an  evaluation  of  their  relative  importance  and  to  allow  diagnostic  experiments 
to  be  designed.  But  the  effects  of  high  gas  temperature  are  not  well  understood  theoretically, 
and  more  data  are  needed  over  a  wider  range  of  parameters.  The  effects  of  convection  are 
important  in  terms  of  the  conflicting  theories,  each  of  which  can  be  supported  by  some  form 
of  experimental  data.  The  effects  of  slender  reentry  vehicle  flow  fields  on  the  solution  of  the 
electron  continuity  equation  need  to  be  examined  more  carefully  in  terms  of  breakdown  criteria 
in  order  to  resolve  the  convection  conflict  and  provide  a  meaningful  definition  of  pulse  break¬ 
down  in  the  presence  of  flow  field  ionization. 
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